Effect of spin excitations on the property of quasiparticles in electron-doepd cuprates 
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It is proposed that the 50 - 70 meV dispersion anomaly (kink) in electron-doped cuprates revealed by re- 
cent angle-resolved photoemission spectroscopy experiments is caused by coupling with the spin fluctuation. 
We elaborate that the kink exists both along nodal and antinodal directions, and both in the superconducting 
and normal state. The renormalized effect for the density of states is also studied and the hump feature out- 
side the superconducting coherent peak is established, consistent with recent scanning tunnelling microscopy 
experiments. 
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Although high-temperature superconductivity in cuprates 
was discovered more than twenty years ago, the mechanism 
of their unusually high critical temperatures has not yet been 
clarified Ol]. An insightful view may be obtained through 
the understanding of the role played by certain collective ex- 
citations by studying the renormalized single-particle prop- 
erties. Experimentally, the Angle-resolved photoemission 
spectroscopy (ARPES) and scanning tunnelling microscopy 
(STM) have been powerful tools for providing the electronic 
structure and probing the interaction of the quasiparticle with 
certain boson modes. 

The superconductivity in cuprates can be achieved by dop- 
ing either holes or electrons into parent antiferromagnetic 
(AF) Mott insulators. One of the most important features in 
hole doped cuprates revealed by the ARPES experiments is 
the slope change of the quasiparticle dispersion (kink) from 
the momentum distribution curve (MDC) |2]. The kink is ob- 
served along both nodal and antinodal directions at the ener- 
gies about 40 ~ 80 meV. In the past few years, the origin of 
the kinks attracted intensive study both theoretically and ex- 
perimentally because it speculated some kind of interaction 
which might act as the mysterious glue for Cooper pairs. Two 
possible bosonic modes, namely, phonon ^ and spin reso- 
nance mode revealed by neutron scattering experiments fl], 
have been proposed to account for the dispersion kink. Theo- 
retically it seems that both electron-phonon interaction jst] and 
the coupling of the spin resonance mode ^ can repro- 
duce the dispersion kink. Unfortunately, this two modes could 
have similar energies, thus it is difficult to distinguish between 
the two. Up to now no consensus has yet been reached. On 
the other hand, STM experiments also identified the existence 
of the bosonic mode in the hole-doped cuprates, but the origin 
is also under debate fl flolfTill . 

In the past few years, more and more attention has been 
turned to the electron-doped cuprates. It is well known that 
the electron-doped materials exhibit different behaviors from 
that of hole-doped ones, namely, they usually have lower su- 
perconducting (SC) transition temperature and narrower SC 
doping range. Therefore, the spin resonance energy is much 
less, namely, only about 10 meV lfT2i[T3ll revealed by the neu- 
tron scattering experiments. On the other hand, the phonon 
energies are expected to be similar to those of holed-doped 
ones. Thus the energies of these two modes are quite differ- 
ent so that their contributions should be easily separated by 



experiments. Moreover, earlier ARPES experiments in the 
electron-doped cuprates did not observe the kink along the 
nodal direction, only the antinodal kink with the energy about 
50 - 70 meV was observed H H Very recently, it 
was reported by several groups that the kinks exist in sev- 
eral families of electron-doped cuprates, both along the nodal 
direction |[l3, H [H, IMl and antinodal direction ^ [Tsll , 
with the energy being 50 - 70 meV. And the kinks depend 
weakly on the doping level and exist even in the normal 
state. Another renormalized effect revealed by the experi- 
ments is the peak-dip-hump structure in the energy distribu- 
tion curve (EDC) iflST . namely, the EDC line-shapes display 
a sharp quasiparticle peak near the Fermi energy Ef along 
the antinodal direction. This peak terminates and is accompa- 
nied with a dip at the energy about 50 meV. The peak width 
decreases when approaching the Fermi energy. A faint hump- 
like feature is also revealed in the nodal direction. Because 
the kink energy is much greater than the resonance energy 
and the spin resonance peak in fact does not exist in the nor- 
mal state. Thus it was proposed that the phonon should ac- 
count for the dispersion kink iTtI [Tsi [T9ll23l . On the other 
hand, a distinct bosonic mode of the energy about 10 meV has 
also been reported by the STM experiment in electron-doped 
cuprates 121]. It is proposed that the mode is caused by spin 
fluctuations rather than phonons. 

In this letter, the spectral function and density of states in 
electron-doped cuprates observed by experiments iflTi [Tsi [T9I 
I20I I21I1 can be reproduced by only considering the coupling 
between the spin excitations and the quasiparticle. We as- 
sume phenomenologically that the spin excitations are from 
spin fluctuations and the retarded Green's function G(k, cS) is 
a function of the bare normal state quasiparticle dispersion ey, 
the SC order-parameter Ak, and the self-e nerg ies S(k, w) due 
to the coupling of spin fluctuations ll22i l23l l24ll . The bare 
normal state quasiparticle dispersion is expressed by, 

fik - -2fi(cos kx + cos ky) - 4f2 cos k^ cos ky 
-2f3(cos 2kjc + cos Ik,,) 
-4f4(cos kx cos 2ky + cos ky cos 2kx) 

-4f5 cos 2kx cos 2ky - to, ( 1 ) 

with ro-5 = -82, 120, -60, 34, 7 and 20 meV. This single- 
particle dispersion was used by Ref. ||25I1 to fit the ARPES 
experiments in electron-doped cuprates flSll . 
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FIG. 1: The imaginary part of tiie spin susceptibility versus the en- 
ergy (X) at the AF momentum Q = (n, n) in the SC state and normal 
state, respectively. 



The SC order parameter is chosen to have t/-wave symme- 
try, namely. 



Ak - Ao(cos ky - cos A;y)/2. 



(2) 



The spectral function of the electrons can be calcu- 
lated from the retarded Green's function as A(k, w) - 
-(l/7r)ImGii(k, oj + iS). Here the Green's function Gij (i, j - 
1,2) is calculated by Dyson's equation in the Nambu repre- 
sentation (2x2 matrix), namely, 

G(k, w + iSy^ = Go(k, a> + iST^ - S(k, w + id). (3) 
The bare Green function of the electron Go is expressed by. 



ito - Sk -Ak 
-Ak ico + Sk 



(4) 



The self-energy due to spin fluctuation is written as |12 

S(k, iu) = ^ ^ ^ 8^X(<1^ itD,„)cr2Go(k - q, ico - /w,„)a=3, 

(5) 

where 0-3 is the Pauli matrix, xiq, i<^m) is the spin susceptibil- 
ity in the random phase approximation (RPA), namely. 



/ro(q, ^) 

1 + Ut^xo{({,CLl)' 



(6) 



Here t/q = f/o(cos + cos ^y) consistent with the t - J type 
model. Xoi^^ w) can be calculated from the Fermionic bubble, 

;kfo(q,w) = ^ Tr[Go(k,/w„,)Go(k-i-q,/w-i-/w„,)]- (7) 



k,/tt>,„ 



In the follo wing pre sented results, we set Uq - 260 meV, 
g - 360 meV |l 27l |28[ I29I1 . The temperatures and gaps in the 
SC and normal states are T - 0.5 meV, Aq = 10 meV and 
T - Tc - 2.3 meV, Aq - 0, respectively. We have checked 
numerically that the main results are not sensitive to the choice 
of the parameters. 

The imaginary parts of the spin susceptibilities as a function 
of the energy o) are plotted in Fig. 1. A sharp resonance peak 
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FIG. 2: (Color online) The intensity plots of the spectral functions 
as functions of the momentum and energy in the SC state (a-b) and 
normal state (c-d), respectively. The left and right panels are along 
(tc, 0) to in, n), and (0, 0) to {n, n) direction, respectively. The solid 
and dotted lines are the MDC dispersions and the bare band disper- 
sions, respectively. 



is seen in the SC state at the energy about = 10 meV. The 
origin of the resonance has been studied intensively ll30l l3lll . 
It arises from a collective spin excitation mode corresponding 
to the real part of the RPA factor (1 + l/qReXi)) equals to zero 
and the imaginary part of the bare spin susceptibility Im^o is 
small. The resonance is absent in the normal state, where the 
peak intensity decreases dramatically and only a low-energy 
broad peak can be seen. While in fact the weight of the spectra 
are at low energies and near AF momentum both in the SC and 
normal states. 

Figs.2(a-d) show the intensity maps of the spectral func- 
tions [A(k, a>)f(a>)] (/(w) is the Fermi distribution function) 
as well as the MDC dispersions in the SC state (up panels) 
and normal state (down panels), respectively. Clear kinks at 
the energy about cat ~ 60 - 70 meV can be seen along the 
antinodal direction. Well defined quasiparticle peaks exist be- 
low the kink energy. At higher energy, the peak intensity is 
small. In the normal state the dispersion kink still exists and 
no qualitative diff'erence can be seen. 

The right panels of Fig. 2 show the nodal data of the spectral 
function. Here, the dispersion kink can also be seen clearly at 
the energy about = 60 meV. The renormalized effect is 
much weaker than that of antinodal direction but does exist, 
which can be seen more clearly by comparing the dispersion 
with that of the bare band. As shown, the renormalized disper- 
sion and the bare one are nearly parallel at high energies while 
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FIG. 3: (Color online) The line-shape A(k, a))/((i») as a function of 
the energy at different momentums in the SC state (a-b) and normal 
state (c-d), respectively. The left and right panels are along {n, 0) to 
{n, 7i), and (0, 0) to in, n) directions, respectively. 
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FIG. 4: (Color online) The real parts of the self-energies Re Z(k, w) 
vs. the energy to at different momentums and the bare band disper- 
sions £k and £ii_Q along the antinodal (a-b) and nodal (c-d) directions, 
respectively. 



the renormalized one bends to low energy at about 60 meV in- 
dicating that the kink is indeed caused by the self-energy. In 
addition, the peak intensity is larger at low energies (w < w^), 
and decreases evidently at the kink energy. Similar with the 
case of antinodal direction, there is also no remarkable dif- 
ference between the spectrum of the SC state with that of the 
normal state along the nodal direction. 

The EDC line shapes along the antinodal and nodal direc- 
tions are plotted in Fig. 3. Along the antinodal direction a 
sharp quasiparticle peak can be seen near the Fermi momen- 
tum Kf following a 50 meV dip, which is consistent with the 
experiment as we mentioned above ifisll . The peak intensity 
decreases dramatically as the momentum is far away from Kf. 
Only a broad peak can be seen at the momentum (;7r, O.IStt). 
Much weaker renormalized effect is obtained along the nodal 
direction, namely, a sharp quasiparticle peak accompanied by 
a high-energy hump-like tail near the Fermi energy. At higher 
energies (w > wt), the peak becomes a little broader while 
it is still well defined and the hump-like tail disappears. The 
curve seems to be symmetric with respect to the peak energy. 
We can also see that there is no qualitative difference of the 
line-shape between the SC state and normal state. 

Our theoretical results reproduce the dispersion kink. 
Though the spin susceptibility shows remarkable difference 
between the SC state and normal state, namely, a sharp res- 
onance peak can be observed only in the SC state, as seen 
in Fig. 1, while the renormalized spectral function show no 
evident difference between the SC state and normal state. Al- 
though here we propose that the spin excitations should be 
responsible for the kink, while in fact the spin resonance phe- 
nomenon is not essential to the kink. In the following we 
demonstrate the origin of the kink and propose that the bare 
band structure and the renormalization by the spin suscepti- 
bility are both important to produce the kink. 

A sound explanation for the dispersion kink can be given 
through analyzing the self-energy. The peak position is deter- 
mined by the pole condition oj - Sk-Re Z(k, w) = in the 
normal state. The real-part of the self-energy is responsible 



for the kink. Performing the summation over ia>„, [Eq.(5)] we 
can rewrite the self-energy as. 



, 1 V r 2t / 1 -/(Sk-q) , 

S(k,w) = — > g^lmx(q,aJi) " T'^'^i' 

TtN ^ J a> - LiJi - Sk-q + 10 

(8) 

where b{a)) is the Bose distribution function. The real part 
of the self-energy is calculated by using the parameters of the 
bare band. The summation over can be written as the inte- 
gral form, l/(7TN)Yjq 47r Jdq An Jt/gk-q/['^ek-q/'^q]- 
The spin susceptibility is peaked at the AF momentum Q and 
very low energy. As a result, approximately, the absolute of 
self-energy | 2 | should have the maxima value and a kink is 
expected near the flat band, namely, VkSk-Q - 0. 

We show the real parts of the self-energies in the normal 
state and the bare band dispersions along the antinodal direc- 
tion in Fig. 4(a) and Fig.4(b). As seen in Fig.4(a), the absolute 
values of the real part of the self-energy | Re 2 | reach the lo- 
cal maximum at the energies about 50 meV and 400 me V. The 
origin of the two peaks can be seen from Fig.4(b), namely, the 
band dispersion is flat (VkSk-Q = 0) at the energies about 50 
meV and 400 meV. As a result, the MDC dispersion has an ob- 
vious kink at the energies about 50 - 60 meV along the antin- 
odal direction. The kink is always there, regardless Im;t'(Q, o)) 
has a resonance peak or not (see Fig. 1). 

The real parts of the self-energies and the bare band disper- 
sions along nodal directions are shown in Figs.4(c) and (d). 
As seen from Fig. 4(d), below the Fermi momentum k < Kf, 
Sk-Q is always greater than zero, so that for negative ener- 
gies the self-energy is quite small. And different from that 
along antinodal direction, there is in fact no obvious peak at 
low energies and the absolute value of Re E is maximum at 
zero energy, indicating that the renormalized effect for peak 
position is prominent at low energies, consistent with the dis- 
persion shown in Fig. 2. As the momentum is away from the 
Fermi momentum, the self-energy tends to be a constant at 
high energies so that the renormalized dispersion is parallel to 
the bare one at high energies, which can be seen in Fig. 2. 
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FIG. 5: The density of state as a function of the energy in the SC 
state. 

For the case of the SC state, because the SC gap is much less 
than the kink energy, and the spin susceptibiHty is still peaked 
at very low energy. Thus in fact the SC gap does not influence 
the kink very much. We have also check numerously for dif- 
ferent gap symmetry, i.e., the nonmonotonic li-wave gap ifisll 
and obtain similar results. 

We now turn to address the renormalized eff'ect of the den- 
sity of states [p(<y) = jA(k, (y)t/k] in the SC state. Fig. 5 
shows the density of states as a function of the energy. The 
SC coherence peaks at the energies +Ao can be seen clearly. 
Outside the gaps the hump-like features exist, and this reveals 
the existence of the bosonic mode, which is sensitive to the 
intensity and energy of the spin resonance peak and located 
at the energies about +(Ao -i- Q.,) with O,- ~ 10 meV, thus the 



spin resonance mode should account for the humps outside the 
gap. This result is consistent with recent STM experiments on 
electron-doped cuprates ll2lll . In hole-doped cuprates similar 
renormalized effect caused by the bosonic mode was also pre- 
dicted theoretically L32 | ] and observed by STM experiments 
very recently ||9l licf We can also see that the renor- 
malized eff'ect in the electron part (negative energy) and hole 
part (positive energy) is asymmetric. The intensity of the SC 
coherent peak is smaller and the renormalized hump caused 
by the spin resonance is also weaker at the negative energy 
part. In fact, the ARPES experiments can only examine the 
electronic structure in the electron part so that the possible 
renormalized eff'ect at the energy Aq + Or for spectral func- 
tion A(k, CO) is hard to detect and also is not obtained by our 
calculation. In fact, very recently the kink at about the energy 
20 me V along nodal direction was reported by Ref . [201 while 
this result was not reported by other groups lll7lll8lll9fl . 

In summary, we study theoretically the effect of the spin 
fluctuation mode on the spectral function and density of states 
in electron-doped cuprates. We have elaborated that the spin 
excitation is able to cause the 50-70 meV dispersion anomaly 
and the hump-like feature of density of states observed by re- 
cent experiments. Thus we present a consistent picture of the 
effect of bosonic mode coupling for ARPES and STM spectra 
in electron-doped cuprates. 
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